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On the basis of a three-dimensional Hubbard model, the superconducting mechanism of Celn3 
under high pressure is investigated by the third-order perturbation theory with respect to the 
on-site Coulomb interaction U. Here we propose the d-wave pairing state induced by antiferro- 
magnetic spin fluctuations. The estimated superconducting transition temperature is lower by 
one order than that in the two-dimensional system for the same value of U/W (W^bandwidth). 
This result is consistent with the difference between Celn3 and CeRhlns. 
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Heavy fermion compound Celna 1 has been confirmed 
to exhibit the properties of unconventional superconduc- 
tivity. For instance, in the 115 In-NQR measurement, 2 no 
coherence peak has been observed in the temperature 
dependence of 1/Ti. The superconducting (SC) state ap- 
pears under the pressure P with a critical value P c = 2.55 
GPa, and the maximum value of the SC-transition tem- 
perature T c max = 0.2 K. The SC state exists near the 
antifcrromagnetic (AF) phase with the ordering vec- 
tor Q — (tt, tt, tt) at Ce atoms. 3 The situation strongly 
suggests that superconductivity should be connected to 
three-dimensional (3D)-AF spin fluctuations. 

In this paper, we assume that the Celn3 system is given 
by the Fermi liquid state at low temperatures and study 
the SC mechanism induced by the wave number depen- 
dence in the effective interaction between quasi parti- 
cles originating from the Coulomb interaction among /- 
electrons. The SC state is considered to be realized on 
the main Fermi surface of Ce 4/-electrons. 4 We pursue 
a possibility of the d-wave pairing state due to AF spin 
fluctuations near Q = (tt,tt,tt). To describe the heavy 
fermions, here a 3D Hubbard model is adopted and the 
effective interaction is evaluated on the basis of the third- 
order perturbation theory (TOPT) in terms of the on- 
site Coulomb interaction 5 . 6 It is concluded that the SC 
mechanism of Celn3 is d-wave pairing induced by 3D-AF 
spin fluctuations near Q=(tt, tt, tt). 

In previous papers, 7 ' 8 two-dimensional (2D) and 3D 
single-band Hubbard models were generally investigated 
within fluctuation exchange approximation (FLEX). Al- 
though they treated different AF spin fluctuations de- 
pending on dimensions, they obtained the same result 
that T c in the 3D system is lower than that in the 2D sys- 
tem. In addition to the characteristics of the AF spin fluc- 
tuation, we treated different wave number dependences 
due to the vertex correction in 2D and 3D dimensions. 9 

In this paper, we focus on Celn3 to discuss the micro- 
scopic mechanism of superconductivity. On the basis of 
TOPT, we investigate the effect of the interaction with 
the complex wave number dependence on the SC mech- 
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anism in 3D Celn3. The calculation based on TOPT in- 
cludes also the normal self energy in this study of Celn3. 
We treat the wave number dependence originating from 
the main Fermi surface of Celn3 and the calculated T c is 
reduced by renormalized Fermi energy Ef- As a result, 
we performed a new analysis of Celn3 in detail on the 
basis of TOPT. Next, we clarify the suppression of T c by 
calculating the dependence of T c TOPT /T c RPA - liko on di- 
mensionality. The change in T c between 2D and 3D sys- 
tems applies to that between 2D CeRhIn 5 and 3D Celn 3 . 
With regard to CeTIn 5 (T=Co, Rh and Ir), Takimoto, 
Hotta and Ueda 10 studied the SC mechanism with a 2D 
orbital-degenerate Hubbard model. They indicated the 
appearance of the <i-wave superconductivity next to the 
AF magnetic phase in the case where the crystal splitting 
energy is large. On the other hand, Nisikawa, Ikeda and 
Yamada 11 investigated the d-wave superconductivity in 
CelxzCoi-aJns with the 2D single-band Hubbard model. 

Here, we apply the single-band model to both Celn 3 
and CeRhlns due to the following reason. The main 
Fermi surface of Celn 3 is a single band as shown by the 
band calculation of Betsuyaku. 4 The 3D single band in- 
duces the wave number dependence of the AF spin fluc- 
tuation, which plays the main role in superconductivity. 
Both Celn3 and CeRhlns have the same kind of Ce com- 
pounds. Therefore, we apply TOPT to the superconduc- 
tivity in the single band of CeRhlns as well as Celn3. 
The main difference between the two materials is the di- 
mensionality of the Fermi surface. We clarify the effect 
of dimensionality on the wave number dependence in the 
analysis based on TOPT, and moreover, we show the 
justification of the quasi-2D model in CeRhlns. 

We explain the formulation in the following. The Hub- 
bard Hamiltonian is given by 

W = -*i 4 a c i+aa + t 2 ^2 4<j c i+b<7 + U ^2 n^nn, 

i,a.tr i.b.cr i 

(i) 

where c\ a is an annihilation operator for a quasi particle 
with spin a at site i, a and i are, respectively, the vectors 
connecting nearest-neighbor and next-nearest-neighbor 
sites in a simple cubic lattice. Transfer integrals t\ and ti 
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Fig. 1. Diagram of the normal self-energy T, n (q,uj I1 ). The solid 
line is the bare Green's function Go- The broken line is the 
Coulomb interaction U. The broken line of U connects only solid 
lines possessing opposite spins. 
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Fig. 2. Diagrams for the effective interaction of singlet pairing 
within the third-order perturbation with respect to U. The solid 
line is the bare Green's function Go- The broken line is the 
Coulomb interaction U. The broken line of U connects only solid 
lines possessing opposite spins. The two external lines have oppo- 
site spins. The effective interaction is divided into the RPA-like 
part and the vertex correction. The latter begins with the third- 
order terms. 



denote nearest-neighbor and next-nearest-neighbor hop- 
ping amplitudes, respectively. U is the on-site Coulomb 
interaction, and ni CT =Cj Ci CT . Here, we consider that the 
parameters t\, £2 and U include a common renormal- 
ization factor 11 for constructing the heavy fermion. As 
shown in Fig. 3, we adjust the dispersion E^ so as to 
reproduce the main Fermi surface of heavy fermions con- 
structed mainly by Ce 4/-electrons in the simple cubic 
structure, 4 leading to 



.Ek = —2<i (cos k x + cos k y + cos k z ) 

+4^2 (cos k x cos k y + cos k y cos k z + cos k z cos k x ). (2) 

Then, we obtain the bare Green's function of the quasi 
particle as Go (k) = 1 / [ioj n — (Ef. — /j,o)] , where k is a short- 
hand notation defined as fc=(k, ui n ), k is momentum and 
uj n = irT(2n + 1) is the fermion Matsubara frequency 
with temperature T. Note that the chemical potential 
Ho for the non-interaction case is determined by the 
electron number n (per site and spin) as n=^2 k Go(k), 
where ^2k = (T/N) En aim N is the number of sites. 
The dressed normal Green's function G(k) is given by 
G(k)=l/[iu n — (E k — /i) — S n (fc)], where £ n (fc) is the 
normal self-energy, given by TOPT with respect to U as 



Sn(fc) = 

^{(7 2 X o(fc-fc')+^ 3 [Xo(fc-fc')+0o(fc+^)]}Go(fc , ) I (3) 

k' 

with 



Xo 



(q) = -J2Go(k)G (q + k), 



(4) 



Here, q denotes a short-hand notation <?=(q, u n ), where 
v n ~2TiTn is the boson Matsubara frequency. Note that 
the chemical potential [i, shifted from hq, is again deter- 
mined by the condition n~^ k G{k). We show the dia- 
gram of the normal self energy in Fig. 1. 

An effective pairing interaction V between quasi par- 
ticles is evaluated using TOPT. Although the origin of 
superconductivity is investigated by total terms in V, in 
order to analyze the role of V in detail, it is convenient 
to divide it into two parts as 



V(k, k') = V RPA (k, k') + Kcrtc X (fc, k'), 



(6) 



where Vrpa represents the terms obtained by the ran- 
dom phase approximation (RPA) and Vvertex indicates 
the third-order vertex correction terms. The RPA-like 
term reflects the nature of simple spin fluctuations, while 
the third-order vertex correction terms originate from the 
electron correlation other than the spin fluctuations. For 
singlet pairing, Vrpa and Vertex are given by 

V R p A (fc, k') = U + U 2 Xo (k -k) + 2U 3 X 2 o(k - k% (7) 

Vvertex(fc, = 



2U 3 ReY,G a (k- 

k" 



-k')(xo(k+k")-Mk+k"))G (k"). 

(8) 

We show the diagrams for singlet pairing in Fig. 2. 

An anomalous self-energy E a is expressed using 
V(k, k') and an anomalous Green's function F(k) 
as Y> a {k)=-Y lk ,V{k 1 k')F{k'). At T=T C , the linearized 
Eliashberg equation including E a and F(k) is reduced 
to the eigenvalue equation, 

(9) 



xi:t(k) = -J2v(k,k')\G(k')\^t(k'). 



When the eigenvalue A becomes unity, the SC state is re- 
alized and T c is obtained. We solve the equation on the 
assumption that E a has singlet or triplet pairing sym- 
metry. We divide the first Brillouin zone into 64x64x64 
momentum meshes and take N{ = 512 for Matsubara 
frequency u n . The bandwidth W (W 3D ~ 12ti) is a nec- 
essary range of uj n for reliable calculations. The range is 
covered with the condition: \W\ < irTNf. To satisfy the 
condition, we calculate in the region with T >0.0037. 

As a result, the dominant symmetries are d x 2_ y 2- and 
<i3 Z 2_ r 2-wave pairings, which are degenerate due to the 
space symmetry of the cubic system. In Fig. 4, we show 
the wave number dependence of the anomalous self- 
energy T, A (q,uj n — irT) for d x 2_ y 2-wave pairing. On the 
other hand, we did not obtain stable solution in the 
present calculations for other pairing symmetries such 
as dry-wave. 

We explain in detail the mechanism of the rf-wave 
pairing, which indicates d x 2_ y 2 or d 3z 2_ r 2-wave pairing. 
To describe the main large-volume Fermi surface 4 (see 
Fig. 3), we choose a parameter set as t% = —0.2 and 
n=0.45, near the half-filling n=0.5. The main Fermi sur- 
face with a nesting property enhances the bare suscep- 
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Fig. 3. Fermi surface in a quarter first Brillouin zone and bare 
susceptibility xofe^n = 0). 



tibility Xoi^i v n), due to the feature of 3D-AF spin fluc- 
tuations near Q=(ir, tt, n), as shown in Fig. 3. The AF 
spin fluctuation originating from the RPA-like term pro- 
vides an advantageous contribution to the eigenvalue A 
for the d-wave pairing, as shown in Fig. 5 . In the case in- 
cluding only the RPA-like terms, A always increase with 
increasing U. However, it is significantly suppressed for 
a large U value when all terms in TOPT are taken into 
account, since the vertex correction terms suppress the 
d-wave superconductivity. 

Next, we show the U dependence of T c in Fig. 6. Here, 
it is emphasized that the evaluated T c in the unit of t\ is 
consistent with the experimental SC-transition temper- 
ature T c =0.2 K in the value T c ~ 0.003. We estimated 
the renormalized T c with the electron effective mass to* 
obtained by the experiment of dHvA. The details of the 
estimation are as the follows. The renormalized band- 
width W is given as zWq with the renormalization factor 
z. As another relation, W is given as W = 12t\ from the 
t2 = —0.2. Here, Wo is a bare bandwidth. The renormal- 
ization factor is defined by z = mo/m* and too is the 
bare electron mass, dispersion £k of eq. (2) at There- 
fore, the relation of W is 12t\ — zWq. The bare band- 
width W is Wq ~ 0.1 Ry=1.58 x I0 4 K in the band 
calculation, z is obtained as z — 1/16 using the cy- 
clotron mass of dHvA. 12 From the relation, the renor- 
malized t\ is obtained as t\ — 0.819 x 10 2 K. By means 
of ti, the calculated value T c ~0.003(ii) at U ^9.0 cor- 
responds to T c =0.246 K, which is near the experimental 
SC-transition temperature T c =0.2 K. Thus, the present 
calculation for T c well explains the possibility of the d- 
wave pairing state in Celn3 . Because the calculated T c on 
the basis of TOPT is near the experimental SC-transition 
temperature, we consider that the system has a strong 
electron correlation U ~ 9.0, which is larger than the 
value of the bandwidth. 

Furthermore, it is quite instructive to consider the 
comparison between 2D and 3D cases. We investigate 
the effect of dimensionality on the constant strength of 
the substantial electron correlation. Here, the strength 
of the substantial electron correlation means the value of 
U in comparison with the bandwidth. Thus, U/W keeps 
a constant value in the change from 2D to 3D systems. 
For the same value of U/W - 3/4 (t 2 = -0.2, n = 0.45, 
W 3D ~ 12*i and W 2D ~ 8ii), we obtain T 2D =0.042 
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Fig. 4. Anomalous self-energy T: A (q,ui n = nT) for d x 2 _ B 2-wave 
pairing. 



in a 2D square lattice and T C 3D =0.0067 in a 3D simple 
cubic lattice. Namely, T c for the 3D system is lower by 
about one order than that for the 2D system. This is 
consistent with the experimental findings that T C 3D «0.2 
K for Celn 3 (3D system) and T 2D «2.1 K for CeRhIn 5 
(quasi-2D system). 13 

We show in Fig. 7 the change in T c t opt/ t rpa between 
2D and 3D systems obtained by the present calculation 
using the dispersion 

Ek = — 2ti(cos k x + cos k y + t z cos k z ) 

+4<:2(cos k x cos k y + t z cos ky cos k z + t z cos k z cos k x ). 

(10) 

Here, the hopping integrals t\ and ti in the c-axis direc- 
tion are multiplied by t z . Namely, the dispersion relations 
with t z =0 and 1 correspond to those for the 2D square 
and 3D simple cubic lattices, respectively. As a result, 
the suppression of T c by the vertex correction is stronger 
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Fig. 5. U dependence of the eigenvalue A obtained by TOPT and 
RPA-like terms. 



0.3 

CD 

~ 0.25 



< 

CL 

o 



0.2 
0.15 

0.1 
0.05 






-pTOPTyjRPA-like 


























0.2 








j-TOPT . 








. 0.15 




PA-li 


T r 




1 0.1 
. 0.05 

° 






























TOPT 














6 o.e 



0.2 



0.4 0.6 

tz 



0.8 



0.016 



TOPT, f 2 =-0.2,n=0.45 



0.014 ■ withi; 
0.012 




Fig. 6. U dependence of T c calculated by TOPT. 



in the 3D system than in the 2D system. By including 
the suppression on the basis of TOPT, the difference in 
T c becomes one order between 2D and 3D systems. In 
addition to this, the decrease in T c is slight with the 
increase in the 3D characteristic in the region near the 
quasi-2D system. Actually, in the range that t z is small 
(t z = ~ 0.4), the change in T c by TOPT is small un- 
der the condition that U /W is constant (see the inset of 
Fig. 7). Therefore, the T c calculated on the basis of the 
2D model is appropriate for CeRhlns, even if the band 
structure of CeRhlns possesses a weak 3D characteristic. 

Here, we explain the different suppression depending 
on dimensionality in the analysis of TOPT. We consider 
the wave number dependence of xo on the k x -k y plane 
perpendicular to the fc z -axis, as shown in Fig. 3. In the 
quasi-2D system with AF spin fluctuations, the same 
prominent peaks exist near Q=(tt, tt) on planes at all 
k z . In the 3D system, xo has prominent peaks around 
Q=(tc,tt) on planes only near k z —ir. On the other hand, 
since the prominent peaks exist on planes at all k z in the 
quasi-2D system, the RPA-like term naturally gives high 
T c in the quasi-2D system rather than in the 3D system. 
This fact has been known in previous studies. 7 ' 8 - 14 On 
the other hand, the complex wave number dependence 
exists except for the region around Q. As opposed to the 
prominent peaks, the momentum space possessing the 



Fig. 7. Change in T C T ° PT /T C RPA between 2D and 3D systems at 
the same value of W/U. The inset shows T c TOPT and T c RPA_likc 
obtained by TOPT and RPA-like calculations. 



complex wave number dependence occupies a larger re- 
gion in the 3D system than in the quasi-2D system. For 
example, in the 3D system, the complex wave number de- 
pendence spreads on planes at k z = and 7r/2 in Fig. 3. 
Therefore, the suppression by the vertex correction in the 
3D system is stronger than that in the quasi-2D system. 
In addition to the characteristics of the AF spin fluctu- 
ation, we obtain T c by including the effect of the com- 
plex wave number dependence in TOPT. The calculated 
T c is near to the realistic SC-transition temperature by 
including each wave number dependence for each dimen- 
sional system. Therefore, we compare the T c obtained by 
TOPT with the experimental SC-transition temperature 
in Celn3 and CeRhhis. 

We summarize our study of the superconductivity in 
Celn3. By means of TOPT based on the 3D Hubbard 
model, we explained the mechanism of the superconduc- 
tivity of Celn3 from a microscopic point of view. It is 
concluded that the SC mechanism of Celn3 is d x i_ y i (or 
d 3z 2_ r 2-wave) pairing, mainly induced by 3D-AF spin 
fluctuations near Q=(tt, it, it). In the present calculation 
including the suppression of T c by the third-order vertex 
correction, T c in the 3D cubic systems is lower by one 
order than that in the 2D square system for the same 
value of U/W, in good agreement with the experimental 
results for 2D and 3D Ce-based heavy fermion supercon- 
ductors, CeRhlns and Celn3. It has been pointed out by 
means of FLEX that the superconductivity induced by 
AF spin fluctuations is suppressed in the 3D system com- 
pared with that in the 2D system 7 , 8 but here we stress 
that the difference is obtained also by means of TOPT in- 
cluding the vertex correction in this paper. Moreover, we 
found that the suppression of T c by the vertex correction 
is stronger in the 3D system than in the 2D system for 
the same value of U /W . A part of this paper is contained 
in the proceeding of ASR2002. 15 The authors thank Dr. 
H. Ikeda, Professo T. Hotta and Professer H. Harima for 
valuable discussions. 
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